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NOTES  ON  THE  APPLICATION  OF  MUNK’S  THEORY  TO  DAMPING -IN-PITCH  OF 
BODIES  OF  REVOLUTION  IN  SUPERSONIC  FLOW 

H.  P.  Hitchcock 

The  Appendix  to  H.  E.  Maloy,  Jr’s  report,  "A  Method  for  Measuring 
Damping-in-Piteh  of  Models  in  Supersonic  Flow,"  has  some  typographical 
errors  that  need  to  be  corrected,  and  also  requires  some  explanation 
to  make  its  meaning  clear. 

In  the  equation  after  (8),  hu"  should  be  "u" .  In  the  next  equation, 

HN  II  WpvqM 

^  should  be  ^  .  In  equation  (10),  the  first  term  of  the  right  member 

should  be  negative  and  the  last  bracket  should  be  placed  at  the  end  of 

the  equation.  The  right  member  of  the  equation  for  5  should  be  negative. 

Incidentally,  there  is  not  distinction  between  this  &  and  the  a  used 

s 

previously;  in  the  main  report,  this  angle  is  denoted  by  a. 

The  cross  sectional  area  S,  introduced  in  equation  (6),  is  a  function 
of  the  axial  distance  x  from  the  center  of  gravity  of  the  body.  At  the 
nose,  x  =  xq  and  S  =  0.  At  the  base,  x  =  x^  and  8=8^. 

The  static  angle  of  attack  is  a  •  But  if  the  axis  of  the  body  has 

s 

an  angular  velocity  of  rotation,  fl,  each  point  on  the  axis  has  an 
additional  velocity  fix,  as  shown  by  (7)»  and  this  produces  an  additional 
angle  of  attack 

Ox 

a  =  *tt—  , 

£0  U 

where  U  is  the  free  stream  velocity.  The  figure  should  show  that  a 

s 

and  0  are  positive  in  a  clockwise  direction. 
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It  is  assumed  that  the  projectile  oscillates  harmonically,  so  that 

2 

a  =  A  cos  tut,  a  =  -  Atu  sin  tut,  a  =  -  Aai  cos  cut. 
s  s  s 

This  assumption  is  convenient,  but  the  same  result  could  be  obtained 

without  it.  Now,  since  fi  =  a  , 

’  s 

x  2  x 

a  =  -Aco  —  sin  a>t,  &  =  -  An  —  cos  tot. 
by  U  '  a)  U 

Now  considering  equation  (10),  the  moment  over  the  whole  body  is 


M  =  pU 


r  &<- 


.2^2 

Anx  sin  cot - ——  cos  tut) 


n 


+  (AUx  cos  cot  -  Ajsc  sin  tut)S*  -  SAtnx  sin  rnt 

bS 


dx. 


where  p  denotes  the  air  density,  and  S'  =  ^  . 

The  moment  of  inertia  of  the  displaced  air  with  respect  to  a  trans¬ 
verse  plane  through  the  center  of  gravity  is  defined  as 


T  =  n 
‘Cg  -  -  J 


/ 


2 

xBt 


x 

n 


The  volume  of  the  body  is 

r *■> 

V  =  I  S  dx. 
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Integrating  by  parts,  we  find  that 
,xh 

xS*  dx  =  x^  -  V, 


r 


\  2  2  f*l. 

xS»dx  =  xbSb-2 


/ 


xS  dx. 


Hence,  the  equation  for  M  assumes  the  form  of  equation  (10). 
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Now  the  spring  constant  is 


k  =  ‘  v)> 


and  the  damping  factor  is 

rr  2 

Making  these  substitutions  and  using  the  expressions  for  a  and  its 

s 

derivatives,  we  obtain  the  equation 

M  =  I  a  +  n  a  +  to. 
eg  s  ro  s  s 


In  the  present  connection,  it  is  appropriate  to  call  p  &  the 

o  s 

damping  moment.  However,  in  the  theory  of  spinning  shell,  the  rolling 
moment  and  the  Magnus  moment  also  tend  to  damp  the  yaw.'1'  Therefore,  it 
is  preferable  to  call  the  term  due  to  yawing  the  yawing  moment. 


Since  Munk’s  theory  neglects  viscosity  and  compressibility,  it 
could  not  be  expected  to  give  reliable  results  at  supersonic  Mach  numbers. 
However,  T.  L.  Smith  has  pointed  out  that  the  theory  which  was  first  used 
for  dirigible  balloons  might  predict  the  right  order  of  magnitude  of  the 
yawing  moment  coefficient  for  small  angles  of  attack,  providing  the 
cross-flow  (normal  to  the  axis  of  the  shell)  is  subsonic.  In  the  present 
case,  it  did,  although  the  theoretical  result  is  less  than  the  experimental 
values,  which  are  given  in  Table  I  of  Maloy's  report . 

A  closer  estimate  is  obtained  by  the  empirical  formula,  which  was 
determined  from  firings  of  projectiles  from  caliber  0.30  to  37mm  at 
velocities  from  2000  to  3050  feet  per  second: 

Kg  =  0.35  L1*5, 


^TCent,  R.  H.  "Notes  on  a  Theory  of  Spinning  Shell"  APG:  BRL  Report  No.  898, 
1954.  See  Eq.  (15).  After  Fowler,  Kent  uses  the  cumbersome  name, 

•yawing  moment  due  to  yawing*. 
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2 

where  L  is  the  length  of  the  projectile  in  calibers.  For  a  5-caliber 
projectile,  this  formula  gives  =  3«91>  which  is  6$  less  than  the  mean 
experimental  value  at  M  =  3-0  and  5$  more  than  the  mean  at  M  =  4.0. 

It  should  be  noted  that  the  experimental  values  decrease  with  Mach 
number.  The  following  table  exhibits  the  averages. 


Mach 

No.  of 

“h 

No. 

Observations 

t; 

-*-•  i 

2.50 

1 

4.37 

3.00 

5 

4.155 

k  nn 

T  •  v  W 

3=72 

4.50 

2 

3.35 

H.  P.  HITCHCOCK 
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Hitchcock,  H.  P.  "Aerodynamic  Data  for  Spinning  Projectiles." 
APG:  BRL  Report  No.  620,  1947-  See  page  13. 
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errors  that  need  to  be  corrected,  and  also  requires  some  explanation 
to  make  its  meaning  clear. 


In  the  equation  after  (8),  "u"  should  be  "U", 


dS" 


In  the  next  equation, 
^  should  be  .  in  equation  (10),  the  first  term  of  the  right  member 

should  be  negative  and  the  last  bracket  should  be  placed  at  the  end  of 


the  equation.  The  right  member  of  the  equation  for  5  should  be  negative. 
Incidentally,  there  is  not  distinction  between  this  5  and  the  ag  used 
previously;  in  the  main  report,  this  angle  is  denoted  by  a-. 

The  cross  sectional  area  S,  introduced  in  equation  (6),  is  a  function 
of  the  axial  distance  x  from  the  center  of  gravity  of  the  body.  At  the 
nose,  x  =  xq  and  S  =  0.  At  the  base,  x  =  x^  and  S  =  S^. 


mil cx  cj+o+-fr»  onrrl  c.  r\-P  of  +  oolr  -fa  rv 
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But  if  the  axis  of  the  hocLy  has 


an  angular  velocity  of  rotation.  Cl,  each  point  on  the  axis  has  an 
additional  velocity  Six,  as  shown  by  (7),  and  this  produces  an  additional 
angle  of  attack 


fix 

=  —  > 


where  U  is  the  free  stream  velocity.  The  figure  should  show  that  a 

B 

and  Cl  are  positive  in  a  clockwise  direction. 
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It  is  assumed  that  the  projectile  oscillates  harmonically,  so  that 

•  ••  2 

d  =  A  eos  o>t ,  Q£  =  -  Ajd  sin  cut,  a  =  -  An  cos  <nt. 
s  s  s 

This  assumption  is  convenient,  but  the  same  result  could  be  obtained 

without  it.  Now,  since  fl  =  &  , 

s 

x  2  x 

a  =  -.Ans  rr  sin  cut,  ds  =  -  An  —  cos  cut. 
a)  U  cu  U 

Now  considering  equation  (10),  the  moment  over  the  whole  body  is 


r  [sc-. 


.22 

M  =  pU  J  |S(-Anx  sin  cot  -  — ——  cos  cot) 


^  I 

+  (AUx  cos  cnt  -  Anx"  sin  cnt)S‘  -  SAnx  sin  cntj  dx, 

where  p  denotes  the  air  density,  and  S'  =  ^  . 

The  moment  of  inertia  of  the  displaced  air  with  respect  to  a  trans¬ 
verse  plane  through  the  center  of  gravity  is  defined  as 


eg 


=  P 


f 


The  volume  of  the  body  is 
V  =  I  S  dx. 


■ 

n 


Integrating  by  parts,  we  find  that 
>Xb 

xS*  dx  =  xbSb  -  V, 


f 


f 


*b 


2  2 
xS*dx  =  xbSb-2 


P 


xS  dx. 


Hence,  the  equation  for  M  assumes  the  form  of  equation  (10). 
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Now  the  spring  constant  is 


k  =  pU^x^  -  V), 


and  the  damping  factor  is 
pQ  =  pUx^Sb. 

Making  these  substitutions  and  using  the  expressions  for  and  its 
derivatives,  we  obtain  the  equation 


M  =  I  a  +  u  a  +  ka  . 
eg  s  'os  s 


In  the  present  connection,  it  is  appropriate  to  call  nQ&s  the 
damping  moment.  However,  in  the  theory  of  spinning  shell,  the  rolling 
moment  and  the  Magnus  moment  also  tend  to  damp  the  yaw.  Therefore,  it 
is  preferable  to  call  the  term  due  to  yawing  the  yawing  moment. 


Since  Munk*s  theory  neglects  viscosity  and  compressibility,  it 
could  not  be  expected  to  give  reliable  results  at  supersonic  Mach  numbers. 
However,  T.  L.  Smith  has  pointed  out  that  the  theory  which  was  first  used 
for  dirigible  balloons  might  predict  the  right  order  of  magnitude  of  the 
yawing  moment  coefficient  Kg  for  small  angles  of  attack,  providing  the 
cross-flow  (normal  to  the  axis  of  the  shell)  is  subsonic.  In  the  present 
case,  it  did,  although  the  theoretical  result  is  less  than  the  experimental 
values,  which  are  given  in  Table  I  of  Maloy1 s  report. 

A  closer  estimate  is  obtained  by  the  empirical  formula,  which  was 
determined  from  firings  of  projectiles  from  caliber  0.50  to  37ram  at 
velocities  from  2000  to  3050  feet  per  second: 

Kg  =  0.35  L1*5, 


Scent,  R.  H.  "Notes  on  a  Theory  of  Spinning  Shell"  APG:  BRL  Report  No.  898, 
1954.  See  Eq..  (15).  After  Fowler,  Kent  uses  the  cumbersome  name, 

•yawing  moment  due  to  yawing'. 
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where  L  is  the  length  of  the  projectile  in  calibers.  For  a  5-caliber 
projectile,  this  formula  gives  =  5-91,  which  is  6$  less  than  the  mean 
experimental  value  at  M  =  5-0  and  5$  more  than  the  mean  at  M  =  4.0. 

It  should  be  noted  that  the  experimental  values  decrease  with  Mach 
number.  The  following  table  exhibits  the  averages. 


Mach 

No.  of 

No. 

Observations 

i  71; 

c  o * 

-*-•  i  ^ 

a. 

2.50 

1 

4.57 

5.00 

5 

4.155 

4.00 

s 

S' 

V7P 
*  1 

l .  r/\ 

*f.?U 

2 

5. 35 

H.  P.  HITCHCOCK 


2  .  ^ 

Hitchcock,  H,  P.  “Aerodynamic  Data  for  Spinning  Projectiles." 
AFG:  BRL  Report  No.  620,  1947-  See  page  15. 
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A  METHOD  FOR  MEASURING  DAMPING- IN- PITCH 
OF  MODELS  IN  SUPERSONIC  FLOW 


ABSTRACT 

A  technique  has  been  developed  to  measure  the  damping- in- pitch  of  bodies 
of  revolution  in  the  Supersonic  Wind  Tunnel  over  a  Mach  number  range  of  1.75 
to  U.5.  The  angular  position  of  the  model  at  any  time  is  recorded  with  an 
oscillograph.  The  aerodynamic  damping  moments  are  computed  from  measurements 
of  the  rate  of  subsidence  of  the  angle  of  attack. 
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I.  INTRODUCTION 


With  the  increasing  demand  to  know  the  damping- in- pitch  of  missiles 
traveling  at  supersonic  velocities  it.  was  .decided  that  instrumentation  should 
he  developed  for  obtaining  the  data  in  wind  tunnels.  Several  requirements 
must  be  placed  upon  the  apparatus  in  order  to  obtain  useful  data.  First, 
the  apparatus  should  be  able  to  measure  small  damping  moments  since  the  models 
are  small  bodies  of  revolution.  Second,  the  tare  damping,  i.e.,  all  damping 
moments  except  aerodynamic  damping,  should  be  as  small  as  possible  and 
repeatable.  Third,  the  model  support  should  not  interfere  with  the  flow  about 
the  model.  Fourth,  the  apparatus  should  fit  within  the  test  region  with  little 
or  no  modification  to  the  tunnel  structure.  Fifth,  it  would  be  desirable  to 
have  oscillating  frequencies  of  the  model  to  be  the  same  as  free  flight  values . 
Sixth,  the  center  of  gravity  of  the  model  should  be  located  at  the  pivot  axis. 
Seventh,  it  would  be  desirable  to  lock  or  clamp  the  model  to  a  fixed  support 
during  start-up  and  shut-down  of  flow  in  the  tunnel.  Additional  requirements 
were  placed  on  the  instrumentation  so  that  the  effect  of  model  center  of 
gravity  location  and  frequency  of  oscillation  on  damping  could  be  obtained. 

The  model  chosen  for  testing  was  the  Army-Navy  Spinner  Rocket,  which  is 
a  spin  stabilized  body  of  revolution.  The  wind  tunnel  model  was  5  calibers 
long  (11.25")  with  a  secant  ogive  nose  2  calibers  long  (4.50"),  a  base  diameter 
of  2.25"  and  a  center  of  gravity  location  of  4.45"  from  the  model  base  (Figure 
l).  It  was  made  of  a  combination  of  aluminum  and  steel  parts  and  had  a  movable 
weight  located  in  the  nose  to  position  the  center  of  gravity  of  the  model  at 
the  prescribed  position. 

This  report  presents  a  description  of  the  apparatus  and  the  technique 
used  to  measure  the  damping- in- pitch  on  a  body  of  revolution  in  a  supersonic 
air  stream. 
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II.  EXPERIMENTAL  PROCEDURE 


A.  Model  Support  and  Recording  System 

Several  methods  of  obtaining  the  damping  data  were  considered  and  the  one 
that  met  the  prescribed  requirements  was  a  system  which  used  a  crossed  flexure 
assembly  to  support  the  model  and  to  cause  It  to  oscillate  about  a  fixed  point. 
The  oscillation  was  started  and  the  damping-  in-pltch  was  then  determined  by 
measuring  the  rate  of  subsidence.  This  system  had  several  distinct  advantages 
which  are: 

(1)  The  tare  damping,  which  is  due  to  internal  friction  of  the  flexures, 
was  repeatable. 

(2)  Strain  gages  mounted  on  thin  flexures  were  able  to  measure  small 
loads  accurately  and  reliably  if  the  stresses  in  the  flexures  and  gages  yere 
kept  below  the  yield  point  of  the  material. 

(3)  The  flexures  provide  the  restoring  moment  required  to  return  the 
model  to  its  initial  position. 

The  system  had  the  disadvantage  that  the  supporting  strut,  which 
was  made  very  stiff  to  insure  that  the  motion  was  essentially  a  simple 
oscillation  about  a  point,  limited  the  length  and  diameter  of  the  model  which 
could  be  used. 

The  original  flexure  assembly  consisted  of  two  sets  of  flexures,  each  set 
containing  two  flexures  mounted  90°  to  each  other  with  one  end  of  the  flexure 
attached  to  the  model  holder  and  the  other  end  attached  to  the  supporting 
strut.  The  thicker  flexures,  set  A  of  Figure  2,  which  were  the  two  outer 
flexures  of  the  assembly,  provided  the  support  and  restoring  force  to  the 
model.  The  thinner  flexures,  set  B,  which  were  located  between  set  A,  were 
used  primarily  as  strain  gage  mounts.  The  ends  of  the  flexures  were  fitted 
into  square  keyways  and  secured  to  their  respective  pieces  with  screws. 

During  the  initial  tests  It  was  noticed  that  after  the  model  oscillated  several 
thousands  of  times  the  damping  rate  had  changed  with  no  change  in  test  con¬ 
ditions.  Upon  checking  the  assembly  it  was  discovered  that  the  flexure  ends 
were  not  fitting  the  keyways  properly  in  that  only  one  surface  of  the  flexure 
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was  making  metal- to-me tal  contact  with  the  keyway.  Thus  the  only  force  to 
restrain  the  flexure  ends  from  moving  was  the  frictional  force  at  the  surface 
of  the  flexure  end  and  the  bottom  of  the  keyway  and  this  force  was  controlled 
by  the  tightness  of  the  mounting  screws.  Since  it  required  checking  the  unit 
quite  often  during  a  test  a  redesign  of  the  flexure  end  connection  was  made 
in  which  the  square  keyway  was  changed  to  a  modified  ”V”  groove.  This  type 
of  end  connection  required  less  checking  but  after  several  tests  were  com¬ 
pleted  the  flexure  ends  would  move  when  the  model  oscillated.  Therefore  a 
new  design  was  developed  in  which  the  flexures  were  made  as  one  unit,  elimi¬ 
nating  the  individual  end  connections.  This  unit,  shown  in  Figure  3,  was 
clamped  to  the  strut  and  model.  In  addition,  the  flexures  were  all  made  the 
same  thickness.  All  the  test  data  in  this  report,  however,  were  taken  using 
the  original  flexure  assembly  with  the  square  keyway  as  shown  in  Figure  4  (the 
new  flexure  unit  was  not  available  for  these  tests). 

Strain  gages  which  were  mounted  on  the  thinner  flexures  wsrs  us8u.  in 
conjunction  with  a  recording  oscillograph  and  recorded  the  angle  of  attack  of 
the  model  with  respect  to  time.  Two  strain  gages  were  mounted  on  each  flexure, 
one  on  each  surface.  Figure  5  shows  the  wiring  diagram  of  the  assembly  used. 

The  oscillation  of  the  model  was  started  by  a  pneumatically  operated 
trigger  mounted  within  the  supporting  strut  (Figure  4).  The  trigger  traveled 
in  an  arc,  engaged  an  adjustable  horn,  rotated  the  model  and  then  swung  past 
the  horn,  allowing  the  model  to  oscillate.  The  horn  was  adjusted  to  allow 
the  model  to  trip  at  a  3,5°  angle  of  attack. 

A  mechanical  locking  device  was  used  to  prevent  overstressing  the  strain 
gages  and  flexures  when  establishing  flow  in  the  test  section  and  to  stop  the 
model  oscillating  if  it  became  unstable .  This  device  was  a  small  cam  attached 
to  the  supporting  strut  and  operated  outside  the  tunnel  with  a  wire  cable. 

When  the  cam  was  rotated  it  engaged  the  model,  locking  it  securely  to  the  strut. 

The  strut  assembly  was  attached  to  the  wind  tunnel  angle  of  attack  system. 
The  locking  cables  and  air  lines  were  brought  out  through  a  side  access  door  of 
the  wind  tunnel.  Figure  6  shows  the  photograph  of  the  assembly  installed  in 
the  wind  tunnel. 
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The  model  was  located.  In  approximately  the  center  of  the  test  section 
and  at  a  mean  angle  of  attack  of  zero  degrees.  It  was  given  an  initial  angle 
of  attack  of  3*5°  and  then  allowed  to  oscillate  about  its  center  of  gravity. 
As  the  model  oscillated  a  trace  showing  angle  of  attack  with  respect  to  time 
was  raadp  and  simultaneously  one- tenth  second  timing  lines  were  recorded  by  a 
recording  oscillograph.  Figure  7  shows  a  typical  amplitude  record  made  of 
one  test. 


The  model  was  tested  at  various  Reynolds  numbers  over  the  Mach  number 
range  of  1.75  to  4.5.  Additional  tests  were  made  at  several  Mach  numbers  in 
which  most  of  the  boundary  layer  of  the  model  was  made  turbulent.  This  was 
accomplished  by  gluing  a  .20  inch  wide  sand  band  on  the  surface  of  the  nose  of 
the  model  located  2.25  inches  from  the  nose  tip. 


B.  Reduction  of  Data  and  Test  Results 

The  motion  of  the  body  may  be  expressed  by  the  linear  differential  equation 

Ia  +  H2a+k2a=0  (l) 

•• 

where  I  a  is  the  moment  due  to  the  angular  acceleration  of  the  body; 

Mg  a  is  the  aerodynamic  and  mechanical  damping  moment;  and  kg  a  is  the 
aerodynamic  and  mechanical  restoring  moment.  Thus  ^  8111(1  4g  are  expressed  as 

^  *  “o  +  “l 

“2  *  k0  +  h 

where  M  and  k„  are  the  aerodynamic  terms  and  M-,  and  k„  are  the  mechanical 
o  v  lx 

ones.  The  solution  for  equation  1  is 


a  =  a„  e-^(cos  to  t  +  -  sin  co  t) 
u  0) 


(2) 


where 


.  h2  2  *2  ,%2 

X  ■  5J  and  <11  -  T  -  <S) 


with  the  initial  conditions  of 


a  =  0  ,  a  =  aQ  when  t  =  0. 
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The  equation  of  the  envelope  curve  or  the  gradually  decreasing  amplitude 
of  this  damped  harmonic  motion  can  be  written 

n2t 

a  =aQe~  ZT  (3) 

Thus  the  total  damping  term  is 


21  . 


a, 


0 


^2  "  At  xoge  a 
n 

At  is  the  time  interval  between  the  two  model  positions  of  aQ  and  an 


00 


The  maximum  amplitudes  ,  ol 


QU  >  •  ■ 

c.  - 


a.  on  the  oscillograph  records 


were  measured  on  a  Telecordex*  machine  and  the  measurements  which  were  in 
digital  form  were  plotted  on  semi-logarithmic  paper.  However,  many  records 
were  measured  with  a  rule  which  had  one  hundred  divisions  per  inch  and  accurate 
data  were  obtained.  Figure  8  shows  the  plot  of  the  record  shown  in  Figure  7* 

A  mean  straight  line  was  drawn  through  the  plotted  points  and  the  slope  of  the 
line  was  used  to  determine  the  decay  rate. 


The  aerodynamic  damping  coefficient  may  be  determined  after  correcting  for 
the  tare  damping.  The  tare  damping  was  measured  by  oscillating  the  model  on 
the  same  support  system  in  the  tunnel  with  no  flow  and  at  various  pressures 
from  76  to  3.5  cm  mercury  absolute  and  recording  the  motion  on  the  oscillograph. 
The  curve  of  the  damping  term  |x^  versus  tunnel  pressure  was  extrapolated  to 
zero  pressure  to  obtain  the  tare  damping.  Since  the  tare  damping  varied  with 
model  oscillatory  frequency,  data  were  obtained  by  changing  the  moment  of 
inertia  of  the  model  and  the  results  are  shown  in  Figure  9*  Also  it  was 
noticed  that  tare  damping  varied  with  amplitude  and  therefore  to  minimize  this 
effect  data  were  obtained  in  the  2.25°  to  1.75°  amplitude  range  for  all  tests 
including  "flow- on"  and  "flow- off". 


The  Telecordex  machine  is  a  two  channel  indicating  and  recording  instrument 
and  is  capable  of  measuring  the  amplitude  of  each  oscillation  and  presenting 
this  measurement  in  recorded  digital  form. 
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With  the  tare  damping  moment  known  the  aerodynamic  damping  moment  may 
he  obtained.  The  equation  is 


M-qG  =  M-2a  ~  ^ia  ' 

The  damping  coefficient  is  defined  as 

^ti  =  ■  (5) 

pd  u 

The  coefficient  versus  Mach  number  for  the  model  is  shown  in  Figure  10. 
Included  on  the  graph  are  the  results  of  the  free  flight  range  data'1'  and  a 
theoretical  estimate  using  the  Zero  Order  Theory  (presented  in  the  appendix). 

Several  tests  were  made  at  Mach  numbers  3-0,  4.0,  and  4.50  in  which  the 
Reynolds  number  was  varied-  The  results  are  presented  in  Table  I.  As  pan  be 
seen  there  was  an  indication  of  increase  of  damping  with  increase  of  Reynolds 
number.  However,  the  total  spread  of  the  data  was  within  the  experimental 
accuracy . 


Additional  tests  were  made  in  which  a  sand  band  was  placed  on  the  nose 
of  the  model  to  cause  the  boundary  layer  to  become  turbulent.  The  difference 

^  w.  A  m  4  —  m  «•«-.  1  -s  4-  *  1-  V  X  «  4-  «  -w  *4  4-  V  -n  m  J  »  *4  ^ U  .  A-  X-  ma  v.  *4  4.  4  ^ m  X. ^  m  *X  4  w.  X  X.  m  a  a 

xii  l  aoco  uc  owccu  one:  ocd  to  wi  on  aiiu  wi  unuu  o  tiaiioi  oxun  uctnu.  xn  onx  cc 

out  of  four  tests  was  less  than  4  per  cent  which  is  well  within  the  experimental 
accuracy  of  the  data  (refer  to  Table  I  and  Figure  10 ).  Figures  11  and  12  are 

oho^OTJfTi"fi,n>i  TxVir'i-f-  ri  c  Tjh-1  r»Vi  IViD  -f- /-\-P  Vinnnrl  d  1  o  vor  ron  -f  Vid  d  +■  +  Vio 

wuwvuv  ngx  uww  nuxwii  n  uuv  wj  v_»J-  wvxiiu.u>x  j  xu.j  v*x  un  uixw  xuwuwx  uw  wuw 

varioxis  test  conditions.  Schlieren  pictures  were  also  taken  to  show  the 
general  flow  over  the  model  (Figure  13). 


The  test  program  was  started  in  tunnel  No.  3  and  considerable  difficulty 
In  measuring  th  damping  rates  in  the  low  Mach  number  range  was  encountered. 
At  Mach  number  1*73  model  oscillations  wane  induced  by  unsteadiness  in  the 
flow,  making  it  impossible  to  obtain  reliable  data.  A  record  of  such  a  test 
is  presented  in  Figure  l4  showing  the  motion  of  the  model  before  and  after 
it  was  tripped.  As  a  comparison  Figure  15  shows  the  motion  of  the  model 


T  C- 

►  f  J 


The  accuracy  of  the  data  depends  upon  the  amplitude  records,  moment  of 
inertia,  and  the  stagnation  pressure  of  the  wind  tunnel.  The  envelopeof  an 
amplitude  record  was  measured  several  times  by  three  individuals  and  the  maxi- 


measured  on  a  torsional  pendulum  spring  electronic  instrument  with  an  accuracy 
of  2  percent.  The  tunnel  stagnation  pressure  was  measured  to  within  _+  2  mm  of 
Mercury  absolute  of  set  pressure. 

C.  Wind  Tunnel 

The  tests  were  conducted  in  the  Ballistic  Research  Laboratories  *  Flexible 
Nozzle  Tunnels  No.  1  and  3*  They  are  continuous  operation,  closed  circuit, 
variable  density  supersonic  wind  tunnels  with  a  Mach  number  range  of  1.25  to 
5.00.  The  test  section  dimensions  of  No.  1  tunnel  are  15  inches  high  and  13 
inches  wide,  for  No.  3  tunnel  are  13  inches  high  and  15  inches  wide  for  Mach 
numbers  of  2.00  to  5*00  and  20  inches  high  and  15  inches  wide  for  Mach  numbers 
of  1.25  to  2.00.  The  top  and  bottom  walls  of  the  test  sections  are  sloped  to 
account  for  the  boundary  layer  growth  on  the  wal 1 s .  The  Mach  number  distri¬ 
bution  along  the  centerline  for  tunnel  No.  1  are  within  +  .008  of  mean  Mach 
number  and  the  flow  inclination  is  within  +  .1°.  For  tunnel  No.  3  the  Mach 
number  distributions  are  within  +  .030  and  the  flow  inclination  is  within  + 
.5°.  A  schlieren  optical  system  and  a  spark  shadowgraph  unit  are  used  to 
make  visual  observations  of  the  flow  about  the  model. 
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III.  CONCLUSIONS 


1.  All  parts  of  the  model  which,  are  to  be  used  as  a  fixed  unit  should  be 
locked  bight  in  order  to  prevent  any  secondary  motion. 

2.  The  mounting  strut  should  be  stiff  and  tight  so  that  the  motion  is 
essentially  a  simple  oscillation  about  the  cross  points  of  the  flexures. 

3.  The  wind  tunnel  flow  must  be  reasonably  steady  so  that  model  oscil¬ 
lations  are  not  Induced  by  the  unsteadiness  in  the  flow. 

k.  The  damping  coefficient  is  a  function  of  Mach  number  and  decreases 
with  increasing  Mach  number  far  the  model  tested. 

3.  The  damping  coefficient  is  not  very  sensitive  to  changes  in  Reynolds 
number  or  whether  the  boundary  layer  of  the  model  is  laminar  or  turbulent. 

6.  Reasonable  comparisons  with  range  data  were  made. 


H.  E.  MALOY,  JR. 
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IV,  APPENDIX 


As  far  as  the  author  knows  no  exact  theory  exists  that  predicts  the 
damping- in-pitch  of  bodies  of  revolution  traveling  at  supersonic  velocities. 
A  theory  referred  to  as  the  Zero  Order  Theory  was  developed  by  M.  Munk  (Bef . 


The  force  per  unit  length  normal  to  the  x-axis  at  any  given  x  is 


dF  d 

(  AllR  \ 

/  r  \ 

IS 

11: 

w 

1°/ 

ft  x  +  u 

sin  a 

s 

(7) 

S  is  cross  Bectlonal  area  of  model 


Since  a  is  small  Eg,.  "  may  be  written 
ft  x 


u 

U 


TT  -+  a 
U  s 


Let  a  =  A  cos  ait 
s 


u  =  <St  =  -A  a>.  sin  ait  and  a 

B  O) 

*  r 
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-  A  o>  rj  sin  cut 
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(  A  cos  cut  -  A  a)  sin  a>t  ~  )  S 


(8) 


d  A  A  dx  A 
dt-cit+3x  dt~3t+U 

Then  equation  8  becomes 
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The  moment  equation  is 


M  =  F  .x 
x  x 


(9) 


da  da 

S  CD 

Substituting  for  ^ —  and  Fx  and  integrating  over  the  entire 


length  of  the  body 


M  =  2  p  U  IxS  A  (D  sin  cut  dx  -  I  A  co  cos  cut 
1  eg 


+  p  U2  A  cos  cut 


*b  % 


S,  -  vol 


-  p  U  cu  A  sin  cut 


V 

'“b 

r 


.  O  I  O  V  /?v 

I  -  C  j  U  A  UA 


u 
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Let  6  =  A  cos  cot 


5  =  cu  A  sin  cut 


C  ...A _ ,i*4- 

U  =  -UU  t\.  CUB  UUO 


Then  equation  (10)  becomes 


eg 


Itpu^s^tpu2 


x,  S,  -  vol 
1)  b 


\  —  / 


(11) 


Consider  the  6  term  which  is  the  damping  moment  of  equation  (ll)  and  let 


i  =  p  II  x  S. 
o  “b  b 


('12) 


Putting  equation  (12)  in  coefficient  form 


2  „ 

*b  Db 


pUd  d 

where  x^  =  distance  from  base  to  eg 


S,  =  cross-section  area  of  base 
b 


Thus  for  the  model  tested 

^  =  3.07 
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Fig.  6  Photograph  of  model  installed  in  tunnel. 
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Fig.  7  Reproduction  of  typical  oscillograph  record. 
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Fig.  10  Graph  of  aerodynamic  damping  coefficient,)^  ve.  Mach  number. 
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Mach  No.*  i.75 
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Fig.l4  Reproduction  of  an  oscillograph  record  showing  flow 
disturbance  in  tunnel  No.  3  at  Mach  1.73  • 


Fig. 15  Reproduction  of  an  oscillograph  record  showing  absence 
of  flow  disturbance  in  tunnel  No.  1  at  Mach  I.75  . 
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